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ABSTRACT 



We examine the wave-functionals describing the collapse of a self-gravitating dust ball 



^ . in an exact quantization of the gravity-dust system. We show that ingoing (collapsing) 

dust shell modes outside the apparent horizon must necessarily be accompanied by outgoing 
^ ■ modes inside the apparent horizon, whose amplitude is suppressed by the square root of the 

csi ; 

00 . Boltzmann factor at the Hawking temperature. Likewise, ingoing modes in the interior must 

cn ' be accompanied by outgoing modes in the exterior, again with an amplitude suppressed by 

Q . the same factor. A suitable superposition of the two solutions is necessary to conserve the 

dust probability flux across the apparent horizon, thus each region contains both ingoing 
and outgoing dust modes. If one restricts oneself to considering only the modes outside 
the apparent horizon then one should think of the apparent horizon as a partial reflector, 
the probability for a shell to reflect being given by the Boltzmann factor at the Hawking 
temperature determined by the mass contained within it. However, if one considers the entire 
wave function, the outgoing wave in the exterior is seen to be the transmission through 
the horizon of the interior outgoing wave that accompanies the collapsing shells. This 
transmission could allow information from the interior to be transferred to the exterior. 
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I. INTRODUCTION 

Hawking's derivation of black hole radiance |l]] was given within a field theoretic context, 
by calculating the Bogliubov transformation of scalar field operators in the Schwarzschild 
spacetime. The physical process responsible for the radiation was, however, generally taken 
to be a quasi-classical tunneling of particles through the horizon. In the tunneling picture, 
pair production occurs near the horizon. A negative energy antiparticle falls in and a 
positive energy particle is radiated out, leading to a decrease in the black hole mass. Thus 
the intuitive picture and the actual computation were not directly comparable until, many 
years later, the tunneling picture was directly used in [Q, ^ to provide a new semi-classical 
derivation of the Hawking effect (see also @-0). Developments and refinements of these 
methods can be found in |f^-[lO|. 



In the tunneling picture one computes the tunneling rate as the exponential of the imag- 
inary part of the classical action for scalar particles in a neighborhood of the horizon. This 
formalism has also been applied to de Sitter space in [|lT] and shown to have connections with 



black hole thermodynamics in [0. Another approach to understanding Hawking radiation 



is via the gravitational anomalies of an effective chiral theory obtained from a dimensionally 
reduced scalar field theory near the horizon of a Schwarzschild black hole [p!3| , |14|. Because 
the computations are performed in a static background geometry, none of these approaches 
can address the processes that occur during gravitational collapse, neither can they take 
into account the effects of quantum gravity. 

In this paper we will examine the behavior of the wave functionals describing the quan- 
tized collapse of a self-gravitating dust cloud with a view to understanding the near horizon 
processes during gravitational collapse. To do so we will use an exact canonical quantization 
of the LeMaitre-Tolman-Bondi (LTB) models [|^ that was developed not long ago by us in 



16Hl8|. The classical geometrodynamic constraints of the dust-gravity system are given in 
terms of a canonical chart consisting of the physical (area) radius of dust shells, the mass 
contained within a shell (the mass function), the dust proper time and their conjugate mo- 
menta. As usual, the system is described by a Hamiltonian constraint and a momentum 
constraint. For the LTB models it is possible to eliminate the momentum conjugate to 
the mass function using the momentum constraint and the resulting simplified Hamiltonian 
constraint is quadratic in the remaining momenta, yielding a Klein-Gordon-like functional 



differential equation, the Wheeler-DeWitt equation, when Dirac's quantization is apphed. 
The quadratic Hamiltonian constraint requires regularization, which was performed on a 



lattice in [|I^. When care is taken to ensure that the momentum constraint is satisfied in 
the continuum limit, the exact wave-functional describing quantum collapse becomes ex- 
pressible in terms of lattice wave functions, each of which describes a collapsing shell and 
is determined by a set of three equations, viz., the Hamilton- Jacobi equation and two addi- 
tional constraints. These make it possible to determine unique, exact solutions everywhere, 
in the interior and in the exterior, as well as the Hilbert-space measure. 

We used these solutions to examine the post-collapse Hawking radiation spectrum by 
computing the Bogoliubov coefficient in the approximation in which a pre-existing black 
hole is surrounded by dust perturbations ||19|, ^]. The mass function was taken to describe 



a central point mass M plus the dust perturbations, which formally mimicked the quantum 
fields in the original Hawking derivation |l|. In this calculation only the exterior wave- 



functional was necessary. Again, the Bekenstein-Hawking entropy [^ ^ of eternal black 
holes could also be determined via a microcanonical ensemble of states [p^-p5| , but for this 
computation of the black hole entropy only the interior wave functional is necessary. 

In the calculations mentioned in the previous paragraph no attempt was made to match 
the exterior and interior wave functionals or the shell wave functions at the apparent horizon 
during collapse. In this paper we will show that matching shell wave functions at the 
apparent horizon leads to an interesting picture, which requires an ingoing wave on one side 
of it to be necessarily accompanied by an outgoing wave on the other side. The relative 
amplitude of the outgoing wave is suppressed by the square root of the Boltzmann factor at 
a "Hawking" temperature that is inversely proportional to the mass contained within the 
shell. ^ Each pair, i.e., an exterior ingoing wave plus an interior outgoing wave or an interior 
ingoing wave plus an exterior outgoing wave is an exact solution to the quantum mechanical 
problem, but neither is physically acceptable because neither conserves the flux of shells 
across the apparent horizon. To achieve shell flux conservation the two solutions must be 



^ Strictly speaking the Hawking temperature, Th = (SttA/)^^, refers to an eternal black hole of mass 
M. Here we show that the temperature appearing in the Boltzmann factor is (AttF)'^, where F is the 
mass function, which is a function of the shell label and is equal to twice the mass contained within it. 
The temperature therefore varies from shell to shell, decreasing as one goes out from the center of the 
collapsing dust ball. We call it the "Hawking" temperature because of its similarity to the traditional 
Hawking temperature. 



suitably superposed so that, on either side of the apparent horizon, the shell wave functions 
become the sum of an ingoing wave plus an outgoing wave whose amplitude is suppressed 
by the square root of the Boltzmann factor at the "Hawking" temperature. 

This can be viewed in one of two complementary ways. If one confines oneself to the 
exterior, ingoing matter appears to have a non- vanishing probability of being reflected at the 
horizon and the ratio of the reflection probability to the absorption probability is precisely 
the Boltzmann factor at the "Hawking" temperature for the shell. However, when the entire 
wave function is taken into account it becomes clear that a different point of view may 
also be adopted. As we will show, the outgoing modes in the interior have a unit relative 
probability of tunneling through the horizon but, because these modes are suppressed by 
the Boltzmann factor relative to the ingoing modes in the exterior, the net effect is that 
the emission probability is still the Boltzmann factor at the Hawking temperature. This 
tunneling from interior to exterior may provide a mechanism by which information from the 
interior is transferred to the exterior. 

Diffeomorphism invariant wave functionals cannot be directly constructed from the shell 
wave functions described above. We discuss how they should be constructed and show how 
they too can be matched at the horizon. Just as the matching of shell wave functions leads 
to a new way to recover the Hawking temperature, matching the wave functionals provides 
a novel way to compute the Bekenstein-Hawking entropy of the end state black hole [^ p2[ . 

We review the essential features of the quantum LTB model and present the exact quan- 
tum states in section II. In section HI we consider the behavior of shell wave functions across 
the horizon. Here we show that an incoming wave in one region must be accompanied by 
an outgoing wave in the other. We also argue that the two solutions must be superposed 
to conserve shell flux across the horizon. We discuss the construction of diffeomorphism 
invariant wave functionals in section IV. We show that the matching of wave-functionals 
implies that the ratio of the transmission probability to the absorption probability depends 
on the mass of the collapsing dust-ball and is given by the black hole Bekenstein-Hawking 
entropy. We close with a brief discussion in section V. 



II. WAVEFUNCTIONALS OF THE QUANTUM LTB MODEL 

A. LeMaitre-Tolman-Bondi Models 

The LTB models describe self-gravitating dust clouds of arbitrary matter distribution. 
The energy momentum tensor is T^,y = e{T, p)UyUy where e{T,p) is the dust energy density, 
U^{t, p) is the dust four velocity, r is the dust proper time and p labels dust shells that form 
the dust cloud. The LTB line element reads 

rfs2 ^ ^^2 _ JM^dp"" - r\t, p)dn^ (1) 

where R{t,p) is the area radius and E{p) is an arbitrary function of p called the "energy 
function". Einstein's equations lead to the expressions 



where F{p) is yet another arbitrary function of the shell label coordinate p called the mass 
function and the negative sign in the equation for drR is required to describe collapse. 

The mass function represents twice the weighted mass contained within the shell of label 
p. If a scaling is chosen so that the physical radius coincides with the shell label coordinates 
p at r = 0, then it can be expressed in terms of the energy density according to 

Fip)= r eiO, p')p"dp'. (3) 

Jo 

The energy function, on the other hand, represents the total initial energy and depends on 

the velocity profile at the "initial" time, v{p) = [drR{T, p)]t-=q, 

2Eip) =v'ip)- - TeiO, p')p"dp' (4) 

P Jo 

In this paper we shall concentrate only on the marginal models, for which E{p) = 0. In this 

case, the solution to @ is 

2/3 



R{t, p)= p 
The central singularity is defined by the curve 



1 3 r^p^ . 



(5) 



f ^ 2p3/2 
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which gives the proper time at which shells achieve zero physical radius. Of course, so 
long as the above equation for r is not satisfied, p = 0, which is also gives i? = by our 
choice of scaling, is not a physical singularity. Various models are obtained from choices of 
F{p)\ even the eternal black hole can be described by the LTB solution in (|T]) by choosing 
F(p) = 2M0(p), where 6 is the unit step function. 

B. Quantization 

The general spherically symmetric Arnowitt-Deser-Misner (ADM) metric, 

ds^ = N^de - L\dr - N'^dtf - R^dn^ (7) 

can be embedded in the spacetime described by ([l|). After a series of canonical transfor- 



mations described in detail in [|16| and [|T8[, this procedure leads to a canonical description 
of the classical black hole in terms of the dust proper time, T{r), the area radius R{r), the 
mass density function r(r) defined by 



F(r) = / rfr'r(r'), (8) 

Jo 

and their conjugate momenta, Prir), Pr{t) and PriT) respectively. The constraints of the 
self-gravitating dust system in these variables can then be given as 

H = p^^+:fpI-^^q 

H, = r'Pr + R!Pr - rP^ ?« (9) 

where J-" = 1 — F/R. 

Dirac's procedure may be employed to quantize this system. The momenta must be 
replaced by functional derivatives with respect to their corresponding configuration variables. 
There are, however, regularization and factor ordering ambiguities that can be encapsulated 
at a formal level by introducing factors of 5(0) into the functional Schroedinger equation 
and writing the Hamiltonian constraint as |jl8| 



Hm[T,RJ] 






^[r,P,r]?«0 (10) 



where A{R, F) and B{R, F) are arbitrary functions of i?, F, which serve to take into ac- 
count the factor ordering ambiguities. The divergent factors 5(0) indicate that the factor 



ordering problem is unsolved and can be dealt with only once a suitable regularization of 
the Wheeler-DeWitt equation has been performed. The second (diffeomorphism) constraint 
in (^) requires no regularization and may be written as 



T , , , + R - 
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nr,R,T]^0 



(11) 



5T{r) 6R{r) y6T ^ 

Progress with the constraints can only be made after a suitable regularization has been 
applied. Below we briefly describe and apply a lattice regularization. 



C. Lattice Regularization 

For solutions of the constraints, we make the ansatz 
^[t,R,T] = ^(°)[F]exp 



^ I drT{r) WiT,R,F) 



(12) 



where W{t, R, F) is some function to be determined. It automatically satisfies the diffeo- 
morphism constraint. The Wheeler-DeWitt equation is second order in time derivatives so 
both positive and negative energy solutions exist, but we will confine our attention to the 
positive energy solutions above. It is worth noting that any functional 

(13) 



vl>[r, R,T] = Ul--J dr T{r) W(r, R, F) 

would satisfy the diffeomorphism constraint provided that W has no explicit dependence on 
the label coordinate r except through the mass function, F{r). We have chosen U = exp 
so that the wave-functional may also be factorizable on a spatial lattice, whose cell size we 
call 0", taking a — t- in the continuum limit. Diffeomorphism invariance requires that the 
continuum wave-functional and all physical results be independent of the cell size. On the 
lattice, the argument of the exponential function becomes |l^, |18| 

(14) 



j dr r(r) W(r, i?, F) ^ a J] T, W(r,, R,, Fj 



where Tj = T{rj), etc. This turns the wave-functional into a product state. 



nT,R,r] = 1[^,{T„R„F,) = H^jfexp 



'-aJ2^,W{T„R„F, 



(15) 



provided that U = exp. 



Before proceeding further it is necessary to define what is meant by a functional derivative 



when functions are defined on a lattice |T^. The defining equations can be understood by 
analogy with the simplest properties of functional derivatives of the functions J{x) 

SJ{y) 



5J{x) 

6 
6J{x) 



5{y-x), 
dyJ{y) = 



and from these definitions follows 



6J{x) 



dyJ{y)4>{y) = ^{x) 



On a lattice we define, for the lattice intervals x,- and x 



]i 



6J{xi 
6J{xj 



Afx,; 






lim -^ 



(16) 



(17) 



(18) 



where rj labels the i lattice site and 5ij is the Kronecker S, equal to zero when the lattice 
sites Xi and Xj are different and one when they are the same. Just as 6{y — x) is only defined 



as an integrand in an integral, so A(xj 






should also be considered defined only as a 



summand in a sum over lattice sites. Hence 

6 



lim 



a 5^ J(r,) = Imi a 5^ 



and 



It follows that 



6J{rj 



a V" J{ri)(t){ri) = lim o" V" A(ri - rj)(t){ri 



b 



6J{xj 



1,. d 

— i- — hm 



a o— i-o dJi 



(19) 



(20) 



(21) 



where Jj = J{xj). This is compatible with the formal (continuum) definition of the func- 
tional derivative. 



D. Collapse Wave Functionals 

When (|T8|) and (0) are applied to the Wheeler-DeWitt equation in (|T0]) and \Ef [r, R, T] is 
taken to be a product state, one obtains an equation describing the wave functions at each 
lattice point W^ 



drf 



+ Ti 



dR] 



+ A, 



d 
dR, 



2r2- 



+ Bi + 






ij. ^ 0, 



(22) 



but there is a further restriction arising from the diffeomorphism constraint. Inserting the 



ansatz in (KB) into (|22|), we find 



a^T] 



dW 



dr^ 









1 



Or! ^ ' dm ' OK 



+ Bj = 0, 



(23) 



which must be satisfied independently of a. This is only possible if the following three 
equations are simultaneously satisfied at each lattice site 






T. 



9W, 



dR 



i. _l_ TT J_ _|_ A . 



1 
■] / •'3 



0, 



3 



B. = 0. 



dR^ 



dR 



■3 J 



(24) 



Moreover, it is straightforward that the Hamiltonian constraint is Hermitean if and only if 



Aj = J^jdRAia{mj\J^j\). 



(25) 



where xxij is the Hilbert space measure. 

Unique solutions to the equations in (plf ) and having the form given in (|15D have been 



obtained in all, even the non-marginally bound, cases ||18[. For the marginally bound models 
the solution for the phase Wj in the exterior, i.e., for shells that lie outside the apparent 
horizon {Rj > Fj), is 



W. 



(±) 



Tj ± 2F, 



,-1 1 
Zi — tanh — 



^jJ 



, z^>l 



(26) 



where Zj = ^RjjFj. The positive sign refers to ingoing waves, traveling toward the horizon 
and the negative sign to outgoing waves, as can be seen from the signature of the phase 
velocity, 

(27) 



T 



2F,z] ' 



keeping in mind that Zj > 1. In the interior, i.e., for shells that lie inside the apparent 
horizon {Rj < Fj), the solution is 



W^> = r, 



}• ± 2Fj [zj — tanh ^ Zj~\ , Zj < 1 
9 



(28) 



but here the the positive sign refers to outgoing waves and the negative sign to ingoing 
waves, travehng toward the central singularity, again as determined by the phase velocity. 
Furthermore, as shown in Appendix B of [|T^, the system in ( p^ determines not only VU,- 
but the Hilbert space measure, m.j, as well. For the marginal models under consideration, 
tTij is regular everywhere and given by 

m, = zj (29) 

upto a constant scaling. 

III. SHELL WAVE FUNCTIONS 

The interior and exterior wave functions must be matched at the apparent horizon, but 
the phases diverge there because the apparent horizon is an essential singularity of (0). We 
can perform the matching, however, if we analytically continue to the complex plane and 
consider the functions 



W^^\t„z,)=< 



Tj ± 2Fj [zj — Arccothzj] , 

(30) 
Tj ± 2Fj [zj — Arctanh^;^] , 



where Arctanh(2) and Arccoth(2) refer respectively to the principal value of the inverse 
hyperbolic tangent and cotangent functions, for which the following identities are well known 

El: 



Arccoth^; = Arctanh- (31) 

z 



and 



Arccothz = Arctanh;2 + < 



f Im{z) < 

(32) 

-f lm{z) > 



The functions 



;(±) 



^onti'Tj, Zj) = Tj ± 2Fj [zj - Arccothzj] , (33) 



which refer to waves in the region outside the apparent horizon, can be transformed into the 
functions 

y^t^ {rj .Zj)= T, ± 2F, \z, - Arctanhz, ] , (34) 
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which refer to waves in the region inside, by employing the property (|32|) above for lra{z) < 
0. First, assuming that the proper time of the dust shell, Tj, remains unaffected by the 
transformation, one finds 



;(+) 



i+)i 



W',:>{T,,z,) = Wi:\T,,z,)-tnF, 



(35) 



and 



Muhr,,z,] 



,(-), 



^L'{rJ,^,)+^^F,. (36) 

Thus an ingoing wave of positive energy in the exterior gets transformed into an outgoing 



wave in the interior according to 



',(+) 



Fij.{+), 



C,^(r„ R,, F,) = e— ^■^^■^|:^(r,, i?,, R 



(37) 



where Uj = crTj/2 is the shell energy. Likewise an ingoing wave of positive energy in the 
interior is transformed into an outgoing wave in the exterior 



',(-) 



F^J.i-), 



<-^(r„ R,, F,) = e---^-^^-^^-^(r,, i?,, P, 



(38) 



One can easily check that the derivatives of the states coincide in the lower half plane as 
well. 

We can now give two separate solutions to the quantum mechanical problem: since the 
exterior, ingoing wave is matched only to an interior, outgoing wave. 



^ 



(1) 



'^■+2^4\/^-t^^'^"V^ 



— ZU)j 



■'i+^^A\/^'^''''^~y^ 



Rj > Fj 



Rj < Fj 



(39) 



will be continuous and differentiable everywhere. Likewise the interior, ingoing wave is 
matched only to an exterior, outgoing wave, therefore 



^ 



(2) 



7-1 LUU-, 

g-^a;,F,g 



'J-2^4\/^-*^'^h~V^ 



'-E^— tanh ^ \ -cT- 
^3 V ^3 , 



Rj > Fj 



Rj < Fj 



(40) 



is also continuous and differentiable everywhere. The problem is that neither of these so- 
lutions is physically acceptable: the wave functions in (|39D represent a flow toward the 
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apparent horizon both in the exterior as well as in the interior, while the wave functions in 
(^) represent a flow away from the apparent horizon, again in both regions. Therefore we 
consider a linear superposition of the two solutions, 



'^■+2^4\/%-t^"h'V^, 



+ AiC^'^^^^'^e 



^j= { 



„ —lUJi 



> " j V J 



+ A^e 



^j-2^4\/^~*^°'^"N/fl7 



''j-2^4\/^"*''°''~'V^ 



Rj > Fj 



Rj < Fj 
(41) 



where Aj are constants to be determined. 



From the Klein-Gordon-like equation (|2^) for shells one obtains a continuity equation. 



dij.J^ = 0, in which we define 



Jj = -«sgn(J^)mjV^*a^^j, Jf = -imjlTjl^pjdR^iljj, 



(42) 



where the components of J^ have been chosen so that the shell current density, J^, is 
continuous across the horizon. We cannot give an interpretation for JJ within this quantum 
mechanical model, but since JJ is time independent it follows that J^ is divergence free, 
therefore constant. In the exterior and in the interior we find 



J. 



R 



2uj{\Aj 



2uj(e- 



-2lTL0iF, 



'^'^ - 1), R. > F, 



lA-l 



, iuj 



Rj < Fj 



(43) 



and so \Ai 



1. Therefore, we take Aj = 1 for every shell. This gives the absorption 



probability, -Pabsj = 1, which is reasonable since there is no barrier for a shell to cross the 



apparent horizon from the exterior. The right hand side of the first expression in (|4l|) is 
the exterior wave function of a shell. Its first term is the ingoing wave that represents 
the collapsing shell. Its second term is an outgoing wave that represents a shell refiection 
at the horizon and an external observer could view the factor i?ext = e~'"'^^^^ in ipout as 
the refiection coefficient |^, ^. From such an observer's perspective an ingoing wave is 
necessarily accompanied by a refiected wave at the apparent horizon and the ratio of the 
refiection probability to the absorption probability is the Boltzmann factor 



Pr 



refj 



abs,j 



^-2nuijFj^ 



(44) 
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at the temperature 

Tp = i2nF,)-\ (45) 

This is twice the "Hawking" temperature. A similar discrepancy between the calculated 
temperature and the Hawking temperature was also noted in early analyses of the canonically 
invariant, null geodesic tunneling formula [^9|-|33|. It was later discovered ||3^^6| that the 



discrepancy arose because an additional contribution, coming from the temporal part, had 
been ignored. We will now show that the same is true in our approach. 

Let us see how an additional contribution from the temporal part comes about in our 
picture by reconsidering our assumption that the shell proper time remains unaffected by the 
analytic continuation. In the marginal models the proper time is related to the momentum, 
Pr, conjugate to the mass density, F, by p 



r = 2Pr ± 2F [z - Arccoth;z] (46) 

in the exterior and 

r = 2Pr ± 2F [z - Arctanhz] (47) 

in the interior, where the positive sign in both cases is for ingoing matter and the negative 
sign for outgoing matter. If Pr (and not r) remains unchanged by the continuation then, 
analytically continuing as before in the lower half plane, we find 

Tout = Tin T i-rtF. (48) 

In the standard computations leading to the Hawking temperature via the Bogoliubov coef- 
ficient, modes are defined with respect to Killing time and then compared at $5+. Although, 
in the midst of the collapsing dust ball there is no time-like Killing vector and therefore 
no Killing time, Pr has the interpretation of one half the Killing time in the Schwarzschild 
spacetime outside collapsing ball and therefore is its natural replacement inside. 
The additional contribution from the proper time in (^S]) implies that 



Wi-)(r„ z,) = Wir)(r„ z,) + 2z7rF (49) 
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so that (|39|) and ( ^OD should read respectively, 



ij 



(1) 



'^■+2^4\/5|-t^'^h-iy^ 



and 



t, 



(2) 









,j+2F,( ./^-tanh"^ ■ '-^ 



^^-2^4\/^-t^'^'^"V^ 



'-ET-— tanh \ -rr- 
^j V J y 



i?j > Fj 



i?, < F, 



i?j > Fj 



Rj < Fj 



(50) 



(51) 



From these we construct the (corrected) shell wave functions by superposition 



V'j 



,+2F, ( ,/^-tanh-i^/^ 



_j_ ^-2ttujjFj^ 



g-27rWjFjg 



-tanh 



+ e 



,,-2F, ( ,/-52--tanh"i ,/-^ 



-tanh 



i?j > Fj 



Rj < Fj 

(52) 

The ratio of the reflection probability to the absorption probability is now the Boltzmann 

factor at the "Hawking" temperature appropriate for the shell. From the point of view of 

the collapse, this "reflected" piece of the external wave function is a purely quantum effect 

necessitated by the existence of an ingoing wave in the interior. This is not the same as 

the tunneling picture of black hole evaporation in which pair production occurring near the 

horizon causes a negative energy antiparticle to fall into the black hole and a positive energy 

particle to tunnel through into the exterior. Here we have considered only positive energy 

solutions. 

The right hand side of the second expression in ( |5^ ) is the interior wave function of a 

shell. Its second term is an ingoing wave that represents a shell continuing its collapse into 

the central singularity. Without recourse to boundary conditions at the center, which would 

necessarily be ad hoc, we see that it must be accompanied by an outgoing wave (the first 

term in the same expression). This outgoing wave comes with amplitude Rint = e~'^'"^^^^ , 

precisely equal to the amplitude for reflection on the horizon. Once again this is a purely 

quantum effect that is required by the existence of an ingoing wave in the exterior. 
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IV. COLLAPSE WAVE FUNCTIONALS 



The manner in which the shell wave functions match across the horizon is essential for the 
construction of diffeomorphism invariant wave functionals describing the collapse as we argue 
in the following section. The superposed wave functions in (|52| ) cannot be used to construct 
diffeomorphism invariant wave functionals describing the collapse, since they are not simple 
exponentials as required by (|13|). Nevertheless, linear superpositions of diffeomorphism 
invariant functionals will be diffeomorphism invariant and we can construct diffeomorphism 
invariant functionals from each of (^) and (^Tj) by taking the continuum limit of the product 
over shells. Thus in the continuum limit (^Dj) turns into 

g-l/dr r w(+'(r,iJ,F) R>F 

^1 = < (53) 



whereas fFTD into 



^, 



(-)f 



g-Tr/drFFg-i/dr r Wi-^(r,R,F) /^ > ^ 



(54) 



.,(-)r 



and because F = F'{r), the amplitude in each case is 



R<F 



drT{r)F{r) = -{F\oo) - F\0)) = 2M^ 



(55) 
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where M is the ADM mass of the collapsing dust-ball. Thus the diffeomorphism invariant 



functional equivalent of the superposed shell wave functions in 



is vl> = vl>i + vl>2 or 



^ 



(-)r 



e-iJdr r W^t>{T,R,F) ^ g-2^Af2g-l/dr T W^->{t,R,F) ^ > ^ 



(56) 



g-2.M2g ^ 



i/dr r Wy'(T,R,F) 



e 2 



JdrT Wi>(T,R,F) Jl^p 



and the ratio of the reflection probability to the probability for absorption is 

-Pref 



-AttM'' 



P. 



abs 



(57) 



where S is the entropy of the end state black hole. A similar matching of the wave func- 



tionals was also found in [|18[. This provides an alternative and direct derivation of the 
Benkenstein-Hawking entropy in the same way as the matching of shell wave functions pro- 
vides a derivation of the Hawking temperature. However, the extent to which this result is 
generic remains to be tested. 
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V. DISCUSSION 

Let us first summarize the main results of tliis paper. We fiave used a canonical quan- 
tization of the marginal LeMaitre-Tolman-Bondi model of a self-gravitating dust ball to 
understand the near horizon behavior of the wave-functionals describing collapse. There 
is generally a factor ordering ambiguity in any attempt to quantize gravity, which can be 
resolved only after a suitable regularization has been chosen. As in previous works, here 
regularization was performed on a lattice. With this regularization exact solutions of the 
Wheeler-DeWitt were obtained. They cannot be directly matched at the horizon, which 
is an essential singularity of the Wheeler-DeWitt equation. Beginning with the shell wave 
functions, we performed an analytic continuation into the complex plane and showed that 
matching each shell's wave function at the horizon requires that an ingoing wave function in 
one region is accompanied by an outgoing wave function in the other. The amplitude of the 
outgoing wave function in each case is given by the square root of the Boltzmann factor at 
the "Hawking" temperature defined by the mass contained within that shell. Furthermore 
the two solutions describing an ingoing wave in one region and an exponentially suppressed 
outgoing wave in the other region must be superposed to conserve the shell flux across the 
horizon. We showed that crossing the horizon involved a rotation of the dust proper time. 
This rotation is required to keep Pr, which is related to the Killing time in the exterior, 
unchanged. 

The probability for an ingoing shell to cross the horizon is unity. There is also an outgoing 
shell in the exterior, which we have called a "reflected" shell but it can also be viewed as a 
tunneling of an outgoing wave in the interior into the exterior with unit relative probability. 
The ratio Prd/ P&hs = e"'^^'^, where u is the infalling shell's energy and Ph is the inverse 
"Hawking" temperature. Matching wave-functionals across the horizon yields the same 
picture, but now Prct/Pahs = ^~^ , where S is the entropy of the black hole that is the end 
state of the collapsing dust. 

The lattice regularization scheme we have used correctly implements diffeomorphism 
invariance in the continuum limit. It also uniquely fixes the factor ordering and it is possible 
to obtain exact and unique solutions to all the constraints. Thus these are the only solutions 
in which states for the dust cloud factorize into (infinitely many) shell states. Other solutions 
would couple the shells, but to find them one would have to find a regularization scheme 
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that is different from the one chosen here. The factorizable solutions are, however, just 
the WKB solutions, i.e., the regularization and resulting factor ordering used here lead to 
states for which the WKB form is exact. If a different regularization scheme, consistent with 
diffeomorphism invariance, can be found then it will produce corrections to the solutions 
given here, in particular to the amplitudes for crossing the horizon and to the black hole 
entropy as computed via (p^). 

Nevertheless the techniques used here should still be applicable. Let us close by illustrat- 
ing the effect of corrections using the results of a recent paper [^ in which the solutions 
have been given up to order l^, where Ip is the Planck length. This will also show off the 
great simplicity of this approach over the traditional Bogoliubov calculation of the Planck 



spectrum. Interpreting the solutions given in [^ as shell wave functions (extended to the 
complex plane), we have 



WW(r,,.,: 



(1 + CII)tj ± 2Fj [zj + ICiy^ - Arccoth^j] , Re{z) > 1 

(58) 
(1 + Cll)Tj ± 2Fj [zj + \Cllz^ - Arctanh^j] , Re(z) < 1 



where C is a constant. Matching now picks up an additional multiple of ^mF coming from 
the multiplier of r,. Thus we find 

Wi^tV., z,) = >vW(r., z,) - 2z7rF(l + ^C/J) 

^^Uir,. z,) = Wl;\r,, z^ + 2z7rF(l + ^C/J) (59) 

This gives the ratio of the reflection probability to the absorption probability as 

-'abs 

which is the Boltzmann factor at the modified "Hawking" temperature 

^^ = 47rF,(l + la^) ^^'^ 

Likewise the corrections to the Bekenstein-Hawking entropy are obtained by matching the 
wave-functionals at the horizon as we have done earlier: 
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In this paper, we have only examined the marginal models for the purpose of illustration. 
The general solutions for non-marginal models have also been worked out both within the 
lattice regularization described in our earlier papers as well as by the authors in [0. We 
intend to examine them within the framework described here in a future publication. 

Acknowledgements 

LCR Wijewardhana was supported in part by the U.S. Department of Energy Grant No. 
DE-FG02-84ER40153. 



[1] S. W. Hawking, Comm. Math. Phys. 43, 199 (1975). 

[2] M.K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000) 

[3] M.K. Parikh Int. J. Mod. Phys. D 13 2351 (2004) 

[4] P. Kraus and F. Wilczek, Nucl. Phys. B 433, 403 (1995). 

[5] K. Srinivasan and T. Padmanabhan, Phys. Rev. D 60 (1999) 24007. 

[6] G. E. Volovik, JETP Lett. 69 (1999) 662. 

[7] M.Arzano, A.J.M.Medved and E.C.Vagenas, JHEP 0509:037 (2005). 

[8] A.J.M.Medved and E.C.Vagenas, Mod. Phys. Lett. A 20, 2449 (2005). 

[9] M. Arzano, Phys. Lett. B 634, 536 (2006). 
[10] R. Banerjee and B. R. Majhi, JHEP 0806:095 (2008); ibid., Phys. Lett. B 662, 62 (2008) 
[11] G.E. Volovik, Int. J. Mod. Phys. D 18 (2009) 1227. 

[12] R. Banerjee, B. R. Majhi, and S. Samanta, Phys. Rev. D 77 (2008) 124035. 
[13] S. Iso, H. Umetsu and F. Wilczek, Phys. Rev. Lett. 96 (2006) 151302; ibid Phys. Rev. D 74 

(2006) 044017. 
[14] K. Murata and J. Soda, Phys. Rev. D 74 (2006) 044018. 

[15] G. Lemaitre, Annales de la Societe Scientifique de Bruxelles A 53, 51 (1933); for an English 
translation see Gen. Rel. Grav. 29, 641 (1997). 
R.C. Tolman, Proc. Natl. Acad. Sci., USA 20 (1934) 410. 
H. Bondi, Mon. Not. Astron. Soc. 107 (1947) 343. 
[16] C. Vaz, L. Witten and T.P. Singh, Phys. Rev. D 63 (2001) 104020. 
[17] C. Vaz, L. Witten and T.P. Singh, Phys. Rev. D 69 (2004) 104029. 

18 



[18] C. Kiefer, J. Mueller-Hill and C. Vaz, Pliys. Rev. D 73 (2006) 044025. 
[19] C. Vaz, C. Kiefer, T.P. Singh and L. Witten, Phys. Rev. D 67 (2003) 024014. 
[20] C. Kiefer, J. Mueller-Hill, T.P. Singh and C. Vaz, Phys.Rev.D 75 (2007) 124010. 
[21] J. D. Bekenstein, Ph.D. thesis, Princeton University (1972). 

J. D. Bekenstein, Lett. Nuovo Cimento 4 (1972) 737; ibid, Phys. Rev. D 7 (1973) 2333. 
[22] J. M. Bardeen, B. Carter, and S. W. Hawking, Comm. Math. Phys. 31 (1973) 161. 
[23] C. Vaz and L. Witten, Phys. Rev. D 64 (2001) 084005 
[24] C. Vaz, S. Gutti, C. Kiefer, T. P. Singh and L.C.R. Wijewardhana, Phys. Rev. D 77 (2008) 

064021 
[25] C. Vaz and L.C.R. Wijewardhana, Phys.Rev.D 79 (2009) 084014. 
[26] M. Abramowitz and LA. Stegun, "Handbook of Mathematical Functions" , Dover Publications, 

Inc., New York, 1972. 
[27] M. Yu. Kuchiev, Phys. Rev. D 69 (2004) 124031; ibid Europhys. Lett. 65 (2004) 445-451. 
[28] M. Yu. Kuchiev, V. V. Flambaum, Phys. Rev. D 70 (2004) 044022. 
[29] B. D. Chowdhury, Pramana 70 (2008) 593. 
[30] P. Mitra, Phys. Lett. B 648 (2007) 240. 

[31] E. T. Akhmedov, V. Akhmedova, D. Singleton, Phys. Lett. B 642 (2006) 124. 
[32] T. Pilling, Phys. Lett. B 660 (2008) 402. 
[33] R-G Cai and L-M Cao, Class. Quant. Grav.26 (2009) 155018. 



[34] T. Nakamura, |[arXiv:0706MT6 



[35] V. Akhmedova, T. Pilling, A. de Gill, D. Singleton, Phys. Lett. B666(2008) 269. 
[36] V. Akhmedova, T. Pilling, A. de Gill, D. Singleton, Phys. Lett. B 673 (2009) 227. 
[37] R. Banerjee, C. Kiefer and B.R. Majhi, "Quantum gravitational correction to the Hawking 
temperature from the Lemaitre-Tolman-Bondi model" arXiv:1005.2264| 



19 



